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The equation of state for five-dimensional hard hyper- 
spheres arising as a weighted average of the Percus-Yevick 
compressibility (|) and virial (|) equations of state is consid- 
ered. This Carnahan-Starling-like equation turns out to be 
extremely accurate, despite the fact that both Percus-Yevick 
equations of state are rather poor. 

PACS numbers: 64.10.+h, 05.70.-a, 05.70.Ce 



Although not present in nature, fluids of hard hy- 
perspheres in high dimensions (d > 4) have attracted 
the attention-ici-a number of researchers over the last 
twenty years.ETo Among these studies, one of the most 
important outcomes was the realization by Fceasier and 
IsbisterEl and, independently—by LeutheusserQ that the 
Percus-Yevick (PY) equationuJ admits an exact solution 
for a system of hard spheres in d = odd dimensions. In 
the special case of a five-dimensional system (d = 5), the 
virial series representation of the compressibility factor 
Z = p/pksT (where p is the pressure, p is the number 
density, ks is the Boltzmann constant, and T is the tem- 
perature) is Z(rj) = J2^=o bn+irf l , where 77 = (n 2 /60)pa 5 
is the volume fraction (<r being the diameter of a sphere) 
and 6„ are (reduced) virial coefficients—The exact values 
of the first four virial coefficients areS't3 b\ = 1, 62 = 16, 
63 = 106, and 64 = 311.18341(2). The fifth virial coef- 
ficient was. estimated by Monte Carlo integration to be 
65 ~ 970.EJ More recent and accurate MonteXIarlo cal- 
culations yield b 5 = 843(4) and b e = 9SS(28)Jid The ex- 
act knowledge of the virial coefficients 61-64 and in some 
cases of the Monte Carlo values for 65 and b§ has been 
exploited to construct several approximate equations of 
state (EOS), several of them being reviewed in Ref. [l4|. 

One of the simplest proposals is Song, Mason, and 
Stratt's (SMS)JZlwhGUD y viewing the Carnahan-Starling 
(CS) EOS for d — 3E3 as arising from a kind of mean- 
field theory, arrived at a generalization for d dimensions 
that makes use of the first three virial coefficients. Baus 
and Colot (BC)Q proposed a rescaled (truncated) virial 
expansion that explicitly accounts for the first four virial 
coefficients. A slightly more sophisticated EOS is the 
rescaled Pade approximant proposed by Maeso et al. 
(MSAV),E3 which reads 



^MSAV 



to) 



1 +P1V+P2V 2 

(1-?7) 5 (1+ 917/)' 



(1) 



to date is the. semi-empirical EOS proposed by Luban 
and Michels.EEl These authors first introduce a function 
£(77) defined by 



Zlm(v) = 1 + 



6 2 r?{l + [6 3 /6 2 -CW64/63] 77} 
1- av)(h/h) v + [Cfo)-l] (64/62)772 



(2) 



Equation (|J) is consistent with the exact first four virial 
coefficients, regardless of the choice of C(v)- The ap- 
proximation £(77) = 1 is equivalent to assuming a Pade 
approximant [2,1] for Z(r\). Instead, Luban and Michcls 
observed that the computer simulation data of Ref. ||, 
{Z s i m (r)i),i — 1,...,8} [cf. Table H|, favor a linear ap- 
proximation for £(77) and, by a least-square fit, they found 
C(?7) = 1.074(16) + 0.350(96)77. Another semi-empirical 
EOS (not included in Ref. |lj) was proposed by Amoros 
et al. (ASV)MI 



Zasv(tj) = J^^+irf + 



n=0 



5770 

Vo-V 



- 1 



(3) 



where p x = (776 - 6 4 )/36, p 2 = (5476 - 116 4 )/36, and 
<7i = (380 — 6 4 )/36. One of the most accurate proposals 



This equation imposes a single pole at the close-packing 
fraction 770 = \^2tt 2 /30. The parameters (3 n = b n — 

5?7q ^ n 1 ' are fixed so as to reproduce the first five virial 
coefficients, while C is determined by a fit to simulation 
dataE By using the presently known values of 64 and 65 
and minimizing £) i=1 [l - Z A sv(r]i)/ z sim{Vi)} 2 one nn ds 
C = 276.88. Finally, PadjLapproximants [2,3] and [3,2] 
have also been considered S3 

All the previous EOS rely upon some extra informa- 
tion, such as known virial coefficients and simulation 
data. On tip, other hand, (approximate) integral equa- 
tion theories! 5 ] provide the correlation functions describ- 
ing the structure of the fluid. From these functions 
one can obtain the EOS, that usually adopts a different 
form depending on the route followed. As said above, 
the PY integral equation has an exact solution for a 
system of hard spheres in odd dimensions. In particu- 
lar, the analytical expressions of the EOS obtained from 
the virial route, Zpy_ v (r]), and from the, . com pressibility 
route, Zpy_ c (ri) are known for d = 5nuSM-i Neverthe- 
less, these two EOS are highly inconsistent with each 
other .El'ta This inconsistency problem is also present with 
lower dimensions (except in the one-dimensional case, 
where the PY theory becomes exact) ^but to a lesser ex- 
tent. This led Freasier and Isbistera to conclude that 
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"the PY approximation for hard cores is an increasingly 
bad approximation as the dimensionality of the system 
grows larger." 

As is well .known, the CS EOS for three-dimensional 
hard sphereal3 plays a prominent role in liquid state 
theory.Ej While originally derived from the observation 
that the numerical values of the known virial coeffi- 
cients came_cemarkably close to fitting a simple algebraic 
expression,^ the CS equation is usually viewed as a suit- 
able linear combination of the compressibility and virial 
EOS resulting from the PY theory,E£l namely 

z cs(v) = a^Z PY _ c { V ) + (1 - a^)Z PY . v (ri) (4) 

with — |. Since, as it happened in the case d = 3, 
both PY routes keep bracketing the true values in the 
case d = 5, til it seems natural to wonder whether the 
simple interpolation formula (Q) works in this case as 
well. This question was addressed by Gonzalez et al.,h3 
who kept the value a^ 5%> = |. The main goal of this Note 
is to propose a different choice for . The virial coef- 
ficients corresponding to Zcs(v) are = a^b^~ c + 
(f — a^)b^ Y ~ v . By using the known values of b^-bg one 
gets, however, conflictive estimates for the mixing param- 
eter c*( 5 ) namely ~ (b 4 - &J Y ~ V )M Y ~ C " ^ Y_v ) - 
0.68, ~ (b 5 - bp- v )/(b^ Y ~ c - 6^ Y " V ) -0.02, and 
~ (6b - &6 Y " V )/(&6 Y " C " 6 6 Y " V ) - 0-40- On the other 
hand, minimization of Yh=i[ 1 ~ Zcsijli) / Zsimi^i)} 2 yields 
a( 5 ) — 0.62. By simplicity, here 1 take the rational num- 
ber = | and propose the corresponding EOS (H). 

Table | compares the MSAV, LM, ASV, and CS EOS 
with available computer simulations a This table com- 
plements Table II of Ref. [14], where .Zmsav> ZasVj and 
Zqs (the latter being proposed in this Note) were not 
included. In general, the accuracy of the EOS with ad- 
justed virial coefficients improves as the degree of com- 
plexity increases. More specifically, the average rela- 
tive deviations from the simulation data are, from worse 
to better, as follows: 4.7% (Z [32] ), 3.7% (Z SU s), 3.0% 
(Zp*), 2.3% (Z BC ), 0.4% (Zmsav), 0.17% (Z LM ), and 
0.15% (Z ASV ). Concerning the two PY EOS, both are 
quite poor, with average relative deviations equal to 7.1% 
(Zpy-v) and 4.2% (ZpY- c )- The most interesting point, 
however, is that the CS-like EOS (0) (with the choice 
a (P) — |) presents an excellent agreement with the sim- 
ulation data (the average relative deviation being 0.3%), 
only slightly inferior to that of the semi-empirical EOS 
i?LM and ^asv- This is especially remarkable if one con- 
siders that only the first three virial coefficients of Zcs 
are exact, a circumstance-also occurring in the case of 
the original CS equation.li!3 The choice = |, til on 
the other hand, yields an average relative deviation of 
0.5%. 

Let me conclude with some speculations. It seems in- 
teresting to conjecture about the existence of possible 
"hidden" regularities in the PY theory for hard hyper- 
spheres explaining the paradox that, although the virial 
and the compressibility EOS strongly deviate from each 



other, a simple linear combination of them might be sur- 
prisingly accurate. Since the adequate value of the mix- 
ing parameter is — § for d = 3 and = | for 
d = 5, it is then tempting to speculate that its general- 
ization to d dimensions might be a^ d ' = (d + l)/2dJlJ so 
that a(°°) = h in the limiLpf high dimensionality, in con- 
trast to other proposals.EJ For d = 7 the above implies 
that, while Zpy. v and Zpy_ c would dramatically differ, 
the recipe (g) with = | could be very close to the 
true EOS. The confirmation or rebuttal of this expecta- 
tion would require the availability of simulation data for 
d = 7, which, to the best of my knowledge, are absent at 
present. 
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TABLE I. Compressibility factor Z as obtained from simu- 
lation and from Eq. (MSAV), Eq. (g) (LM), Eq. (|) (ASV), 
and Eq. (§) with a ( ^> = § (CS, this work) 



pa 5 


Simulation a 


MSAV 


LM 


ASV 


CS 


0.2 


1.653(1) 


1.653 


1.653 


1.653 


1.653 


0.4 


2.624(3) 


2.617 


2.618 


2.618 


2.616 


0.6 


4.008(6) 


4.003 


4.009 


4.007 


4.000 


0.8 


5.997(11) 


5.964 


5.986 


5.979 


5.964 


1.0 


8.748(16) 


8.720 


8.770 


8.758 


8.731 


1.1 


10.523(20) 


10.488 


10.553 


10.548 


10.510 


1.15 


11.589(22) 


11.490 


11.560 


11.561 


11.520 


1.18 


12.217(24) 


12.133 


12.204 


12.219 


12.168 


a Ref. IBI 
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